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Abstract 

A strong electron-phonon interaction in a metal increases the elec- 
tron density of states in the vicinity of the Fermi energy dramatically. 
This phenomenon is called electron-phonon mass enhancement. In 
this paper the question is investigated whether the mass enhancement 
can be manipulated in multi-layers of two metals with strong and 
weak electron-phonon interaction. A rich behavior is observed for dif- 
ferent thickness ranges of the layers. For thin layers one observes a 
rather homogeneous averaged enhancement. However, for an interme- 
diate thickness range the mass enhancement is highly anisotropic, i.e. 
direction dependent, as well as position dependent. For large layer 
thicknesses one obtains the bulk behavior for each metal. 

PACS: 63.20.Kr, 73.21.-b, 74.45.+C 

1 Introduction 

A strong electron-phonon interaction (EPI) alters the electronic properties of 
a metal rather dramatically [3| , jT] . It enhances the electron density of states 
in the vicinity of the Fermi energy. The overall strength of the EPI is summed 
up in a parameter A. Lead is a good example for a metal with strong EPI 
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having a A value of A « 1.6. As a consequence the electron density of states 
at the Fermi energy is enhanced by the factor Z — (1 + A) « 2.6. In addition 
the Fermi velocity v* F = Kkp/fn* is reduced by the factor Z corresponding to 
an enhancement Z of the mass. (Therefore this effect is often called "mass 
enhancement"). Since metals with strong EPI are generally superconducting 
a number of superconducting properties are also enhanced, for example the 
upper critical field B c2 . 

Although there have been many experiments which investigated the prop- 
erties of double and multi-layers of a metal S with strong EPI and a metal 
W with weak EPI, the author is not aware of any theoretical investigation 
of how the contact between S and W influences the mass enhancement in S 
or W. Such an investigation is the goal of this paper. 

In a metal with a strong electron-phonon interaction the Fermi surface is 
not sharp even at zero temperature. This is shown in Fig.l. The occupation 
below the Fermi momentum (we discuss here the case of clean metals) is less 
than one and above the Fermi energy the occupation is not zero but finite. 
At the Fermi energy the occupation drops by the value Z" 1 . The reason for 
this distribution is the following. An electron k' below the Fermi energy can 
virtually emit a phonon (q, A) and make a transition into a satellite state 
k" above the Fermi energy. This process does not fulfill energy conservation 
since the satellite state c£„a* A has an excess energy of AE = e k " + huj q ^ — Sk'- 
The satellite causes a finite electron occupation of the states k" above the 
Fermi surface and a finite hole occupation of the states k' below the Fermi 
surface. 

In Fig.l an electron is introduced into a state k directly above the Fermi 
energy. It changes the occupation of the k -state only by Z^ 1 < 1 because 
this state was already partially occupied, and after the introduction into the 
state k part of the electron makes virtual transitions into the other states 
above the Fermi energy. This reduces the occupation of the state ko. Part 
of the electron is smeared over an energy range of Jtwd in combination with 
single virtual phonons. 

Therefore mass enhancement in an electron state k above the Fermi 
energy has two contributions: 

• The pre-occupation of the state k is larger than 0: The state k is 
already partially occupied before an additional electron is introduced. 

• The post-occupation of the state k is less than 1: After the intro- 
duction of an additional electron the state is not completely occupied 



2 



because the electron makes virtual transitions into the other states 
above the Fermi energy. 

Both contributions are roughly equal. 




Fig.l: The occupation as a function of momentum (energy) for 
a strong-coupling metal. At the Fermi momentum (energy) the 
jump in occupation at zero temperature is 1/Z. 

Since the occupation of the state ko changes only by Z~ x the energy (of 
the quasi-particle in k ) is only £k /Z. Therefore the quasi-particle energies 
are closer together and their density of states is enhanced by the factor Z. 

In appendix Al a short review is given of the treatment of the electron- 
phonon mass enhancement with Green functions. The discussed physical 
interpretation of the mass enhancement in terms of pre- and post-occupation 
is sketched. However, the Green function method cannot be easily expanded 
to multi-layers. Therefore in this paper I use self-consistant perturbation 
theory. 

In part II the electron-phonon interaction and matrix element in multi- 
layers is derived. In part III the amplitudes of the electron-phonon satellites 
are discussed. In a multi-layer these satellites interfere in real space. This is 
discussed in part IV. In appendix A2 the ground-state wave function in the 
presence of electron-phonon interaction is derived in self-consistent perturba- 
tion theory. Finally appendix A3 treats the ground state plus an additional 
electron above the Fermi energy. 
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Even in a homogeneous strong-coupling metal the derivation of the re- 
lation Z — 1 + A is a complex and extended calculation. The multi-layer 
the system is highly inhomogeneous, and I attempt a more modest goal. I 
calculate the occupation of the electron states in real space in self-consistent 
(Brillouin-Wigner) perturbation theory. This yields an expression for the po- 
sition and momentum dependent mass enhancement. A numerical evaluation 
is left to future work. 

Abbreviations: EPI = electron-phonon interaction, EPME = electron- 
phonon matrix element. 

2 Electron-Phonon Interaction in Mult i- Layers 

Now the question is: What happens in a double layer? Does the finite 
occupation above the Fermi energy leak from the strong-coupling metal into 
the weak-coupling metal? If this is the case then one finds an enhanced 
density of states also in the normal metal (because an additional electron 
above the Fermi energy does not change the occupation by one). To analyze 
this question we investigate a very simple model of a double layer. 

2.1 Model 

• We consider a double layer of a strong-coupling metal S and a weak- 
coupling or normal metal W. To keep the analysis simple a primitive 
cubic lattice is used for both films with the same lattice constant a. The 
films are parallel to the x-y-plane. The metal S consists of N z layers 
and lies in the range < z < d s = N z a, and the metal W has Nf layers 
in the regime d s < z < d s + d w = L z = N z a, where d w = N™a and 
N z is the total number of layers in the z-direction. For normalization 
reasons the double film has a finite but large extension in the x- and 
y-direction with the lengths L x = N x a and L y = N y a. Furthermore 
we use periodic boundary conditions in all three dimensions. This 
represents a multi-layer of S and W. (We expect that this is equivalent 
to an isolated double layer with half the thicknesses d s /2 and d w /2). 
The total number of atoms is iV = N x N y (N s z + N?). 

• The electron density in both films is the same and the electrons behave 
as free electrons. The electron states are given by the wave number 
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k = (k x , k y , k z ) with the quantization ki = i^27r/Lj, (i = x, y, z). (The 
quantization in the z-direction extends over the total thickness L z = 
d s + d w ). 

• Furthermore we assume that the elastic properties of both metals are 
identical. Therefore the phonons propagate through both films without 
scattering. The same quantization as for electrons applies to the wave 
vector q = (q x ,q y ,qz) of the lattice oscillations q,i = fii2n/Li, (i = 
x,y, z). Here /ij lies in the range —N i /2<(j Ji < iVj/2. 

Next we derive the electron-phonon matrix element (EPME) <?k 2 -ki,q,A f° r 
the transition of an electron from a state kx into k 2 by absorbing a phonon 
(q, A) or emitting a phonon (— q, A). The EPME is derived in a number of 
textbooks pQ, [2]. Therefore only the essential results are given here. 

The atoms oscillate with the amplitude u n , and their position is given by 

r n = R n + u n 

where R n = na is the average position of the atom n with n= (n x ,n y ,n z ) 
being integer. 

In the lower strong-coupling film for < z < d s the atoms have an 
(somewhat artificial) electron potential 

U(t) = V(t-t u )-V(t-R u ) 

In the upper weak-coupling film for d s < z < d s + d w the atoms have zero 
electron potential. For zero displacement of the atoms the electron potential 
vanishes in both films. Therefore we treat the electrons in both films as free. 
The potential due to the lattice oscillations is 



U (r) « E ns K • V) V (r - R n ) = i£» s E P K ' P) V p e^-^ 

where the sum over n<, represents the sum over all atoms in S. The displace- 
ment u n of the atom at R n is expressed in terms of the annihilation and 
creation operators Qq' t \, a*, ^ for the phonons (q', A) and (— q', A). 
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where e q / A is the unit vector of polarization A. This yields the electron- 
phonon interaction hamiltonian 

H e-p = Ski,k 2 ,q',Afi , k 2 -k 1 ,q',ACk 2 Ck 1 (o q ',A + a *-q>,\) (!) 

with the electron-phonon matrix element g^-k^q'.A 

ffka-id.q'.A = ^v/iV ( 2Ma; 1 ((k 2 -ki) • e q / )A ) V k2 _ kl S f (q'- (k 2 -ki)) 

(2) 

where the structure factor 5/ (Q) [with Q = q'— (k 2 — ki)] is given by 

s f (Q) = ^E n / QR " = ^E^g/^a^aES 1 ^- 0) 

1 1 — e l Q^ aN i 

= Y.G„,G y 5 Q*,GjQ v ,Gy-pj- 1 _ e iQ z a 

where G is a reciprocal lattice vector. 

This means that there is no conservation of the z-component of the lattice 
momentum for the electron-phonon processes. For comparison we denote the 
electron-phonon matrix element for the pure metal S as fl£ 2 _ kl q / \- Its struc- 
ture factor (Q) = Eg^Q.g w hh Q = q'— (k 2 — k x ) fulfills conservation of 
lattice momentum 

^k 2 - kl ,q',A = iy/N \ 2 Mu} ) (( k 2" k l) ' e Q',A) ^-k^G^-k^q'+G (4) 

For those electron-phonon processes in the multi-layer which conserve the 
lattice momentum, i.e., when k 2 — ki= q'+G, one obtains 

N s z o d s 

0k 2 -ki,q',A - ^k 2 -ki,q',A ~ ^ + ^ 0k 2 -ki,q',A 

The weight of the EPME which conserves lattice momentum is reduced by 
the factor d s j (d s + d w ) . 



3 Electron-phonon satellites 

At T = in the absence of the electron-phonon interaction all states within 
the Fermi sphere with k < kp {kp is the Fermi wave number) are occupied 
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and all other states are empty. We denote this state as |^ ) — II c t> l^o) 

k'<k F 

where |$o) is the vacuum. We select from the full Fermi sphere an occupied 
state k'. (In the following the states such k will be sometimes denoted by 
their creation operators such as c£) . Furthermore we choose from the phonon 
spectrum a phonon state (q, A). The state k' can make a transition into a 
state k" (above the Fermi surface) and create a phonon (q, A). Such a state 
can be described as a k"-electron - k'-hole plus one phonon (q, A). We denote 

its amplitude as ak",k', q - The resulting ground state |\I/o) i n the presence of 
EPI is derived in the appendix A2 in self-consistent perturbation theory. One 
obtains for the amplitude of the satellites 

ak "' k ' q ' A " fo£,-^-M (5) 

where T]y is given in appendix A2 by the self-consistent equation JT2J). 

This means that a state c£, with the occupation "1" generates satellites 
c k" a q,A with the (relative) occupation |ai<",k',q,A| 2 - After normalizing the total 
electron state one obtains for the occupation of an electron state kg below 
the Fermi energy and an electron k ' above the Fermi energy 

n (k' ) = — (6) 

1 + (Z)k",q,A \ a k"K v d,\\ 

n (k ') = 



Ek',q,A 


ak' ',k',q,A| 


2 


1 + (Ek',q,A I 


«kJ,',k',q,A 


1 



The step at the Fermi energy is reduced, because (i) a state kg above 
the Fermi energy is partially occupied and (ii) a state k(, below the Fermi 
energy is partially empty. As we discussed above both effects contribute to 
the electron mass enhancement. 



4 The interference of the satellite wave func- 
tions in a multi-layer 

In a multi-layer one expects that the occupations n (k! ) and n (kg) in the 
ground state for energies below and above the Fermi energy are position 
dependent. For simplicity we ignore umklapp processes with G / 0. Then, 
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in a homogeneous metal, one obtains for a given k'-state within the Fermi 
sphere and a given phonon state (q,A) just one satellite state c£„a qA with 
k" = k'— q. This is different for the double layer. Here k" can take many 
possible values, 

k =k-q+v z = k + vg 

where kf( = k' — q and g = , 2 \ z. We call the states c£„ , a* * with the 

amplitude p_ q+I/g! _ q ^/ (r/ k / — £kj;+z/g — ^) a family of satellite states. Each 
member of the family has the same electron hole k'and the same phonon 
(q, A). The members only differ in the wave number of the state k" by z/g. 
The states of a family are coherent and can interfere. Their phases differ by 

exp(ivgz). This yields a modulation of the amplitude of 

real space. 

The total amplitude of the family 

states in real space is 



c k» c k'« q ,A*o > in 



c kJ| Ck ' a q,A^o ) of hole-electron-phonon 



[Vk.' ~ £ k'J+vg ~ nuJ ) 



With Q z = —q z — (k" — k' z ) = —gu we obtain with equ.(j2} 

1 /2 

- iVN ( — — ) V ((~q + y g) " e q^) ^-g+^g 1 1 - exp (-ivgN s z a 

(7) 

Before we evaluate this result in more detail we consider two extreme 



e wgz 



cases: 



1. Only the term with v = contributes. Then we have 

^ * q \ (r) = ^ f *) (-q-^)^jg 
* ' q ' V 2iVH, A ; (tj* - £ k » - tuv) N z 

_ n° g y q , A 

where p° _ q A is the electron-phonon matrix element for the homoge- 
neous metal S. This wave function has constant density in the metals 
S and W. The electron-phonon matrix element of the metal S is aver- 
aged over the total double-layer thickness and reduced by the factor of 
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d 3 / (d s + d w ). The satellite state is equally distributed over both films 
and partially blocks the state k[J in both films equally. This means that 
both films have an identical enhancement factor. 

2. If the contribution of the term z/g can be neglected in the expressions 
{Vk' ~ £ k>J+v S - hu) and ((-q + z/g) -e qjA ) V_ q+!/g then we obtain 



-, i 2-KivNl i 

„,„ i i-^>(-^r) (2™ 

* 2^v=-NJ2+1 AT / o \ eX P T Z 



) 



This yields 



1 1-exp ^ / 27ril/ 



-2 



where S (z) is essentially a step function which is equal to one in the 
strong-coupling metal S and zero in the normal metal W. In Fig. 2 the 
function S (z/a) is shown for a multi-layer with iV* = 10 and N™ = 8. If 
the conditions for case (2) are fulfilled the state k" has its full amplitude 
in the metal S and its amplitude is essentially zero in the normal metal 
W. In this case we expect the full mass enhancement in S and no mass 
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enhancement in W. 




Fig. 2: The amplitude of the electronic part of the 
electron-phonon satellite (with the wave function e~* k o' r ) 
is restricted to the strong-coupling metal S and vanishes 
in the normal metal W. 

For a realistic evaluation we consider a sandwich of 10 atomic layers of 
S and 8 atomic layers of W. For S and W we use the electronic density 
and Debye temperature of Pb with e F — 9.5eV, kp = 1.6 x lO 10 ??!" 1 and 
©d = 90K but assume a simple cubic lattice with a = 3.28 x 10 _10 m. The 
vector g has the value of 1. 9 x 10 9 m _1 . 

Fig. 3 shows a typical electron-phonon process. The electron in the state 
k' below the Fermi surface emits a phonon (q, A) and makes a transition 
into a state k" above the Fermi surface. The amplitude in the satellite 
state is proportional to the inverse energy denominator (77^ — e^r — Huj) 1 . 
Therefore the main contribution is from the regime were |£k"| , l^k'l << h^D 
and the states k' and k" lie close the Fermi energy. 
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Fig.3: The /c z -planes of a double layer with periodic boundary 
conditions. Both metals have the same electron band structure. 
The films have strong and zero electron-phonon interaction. 
The arrows show transitions from state k' to k" with a virtual 
phonon (q, A). For given k', q there are several final states 
k" = k'— q+z/27rz/ (d s + d w ) permitted. The amplitude of these 
states are coherent and interfere. 

In one scenario the main satellite state k' ' lies close to the z-direction. In 
this case the sub-states k" = kg + fg, with negative v are occupied and not 
available. The sub-states with positive v lie above the Fermi energy by an 
energy of 5E (z/) = ve F -^- « 0.06z/e F . Since the Fermi energy corresponds to 
a temperature of about 1.1 x WfiK the next sub-state k" = kg + g lies above 
the Fermi level by an energy corresponding to 6600if . This is very large 
compared to the Debye temperature of 90fC Because of the large energy 
denominator [r]^ — s^" — fkj) the sub-states with k" = kg + z/g can be 
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neglected. In this direction only the main state kg contributes, i.e., only the 
value is — contributes as discussed in case (1). It is remarkable that one can 
increase the film thicknesses by a factor of 100 before the energy SE (is = 1) 
is of the order of the Debye energy. 

The situation is different when kg lies in the x-direction. This is shown 
in Fig. 3 on the left side where k ' point in the negative x-direction. Here the 
states k" = k ' + isg are available. Their energy separation from the state k ' 
is given by 5E (is) = (hug) 2 / (2m) = (vg/k F fe F ss 0.0Uu 2 e F « u 2 15A0K. 
Again this value lies considerably above the Debye temperature of 90K. 
However when we increase the thicknesses d s and d w by a factor 10 then 8E 
reduces to IbK and the sub-satellites have to be included in the calculation. 

This background occupation (at T = 0) yields about one half of the mass 
enhancement. From the dependence of the energy separation SE (is) on is 
and g we obtain the following results as a function of the total thickness 
(d s + d w ) = N z a: 

• N z < 50: The electron-phonon matrix element is reduced by the factor 
p — N s z / (N s z + Nf). Both films are equally enhanced but the enhance- 
ment factor is reduced by p 2 . 

• 50 < N z < 300: The occupation of a state k'J depends critically on the 
direction of k '; for k ' parallel to the z-direction its occupation is the 
same in S and W. Here one still has an averaged EPME. For k ' parallel 
to the film plane the occupation of k ' in S takes its full value while it 
approaches zero in W. The enhanced density of states in both films is 
highly anisotropic. 

• N z > 300: The occupation of the state kg approaches the individual 
value for the two metals S and W. 

The second half of the mass enhancement is due to virtual electron- 
phonon processes which start from the quasi-particle state k with final states 
k", (q, A) where the sum goes over all free electron states k" above the Fermi 
energy. As before the weight of these processes in S and W depends strongly 
on the final state. However, since now we have to sum over the final states 
k" the resulting anisotropy for the state k is strongly reduced. One obtains 
as before the thin film regime for N z < 50 where the EPME is averaged over 
the two metals and the bulk limit for N a > 300 where the individual bulk 
mass enhancement are reinstated. In the intermediate thickness range one 
has a slow transition between the two extremes. 
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5 Conclusions 



In this paper mass enhancement is investigated in double and multi-layers 
of two metals with strong and weak-coupling electron-phonon interaction. 
The mass enhancement is due to the fact that an electron injected into a 
state ko above the Fermi energy changes the occupation of this state by 
less than one. This is for two reasons: (i) the state ko was (even at T = 
0) already partially occupied due to electron transitions from the occupied 
Fermi sea into the state k emitting a virtual phonon (q, A) and (ii) the 
injected electron makes transitions into states k" above the Fermi surface 
emitting virtual phonons (q, A) and reducing the occupation of the state 
k . In multi-layers one has modified electron-phonon matrix elements and a 
finite quantization of electron and phonon states perpendicular to the film 
planes. As a consequence one has to consider interference between electron- 
phonon processes which start from the same initial electron state k' (k' < 
hp), emit the same phonon (q, A) and yield a superposition of different c£„ 
for the final electron state. After the emission of the virtual phonon the 
electron no longer has a well defined momentum. This interference yields 
a spatial dependence of the pre-occupation of the electron state k . This 
results in very interesting properties of the mass enhancement in multi-layers 
of metals with strong and weak-coupling electron-phonon interaction. In the 
thin-film limit the electron-phonon matrix element is averaged over both 
films. In an intermediate thickness range 50% of the mass enhancement 
in each film depends strongly on the direction of the electron momentum 
k . In both films the mass enhancement approaches the bulk value in the 
direction parallel to the film planes while perpendicular to the films one 
obtains an averaged mass enhancement. This will cause a rather anisotropic 
propagation of the conduction electrons parallel and perpendicular to the 
films. In strong-coupling superconductors in contact with normal films it will 
influence the boundary condition between the films and as a consequence the 
superconducting transition temperature of the double or multi layer as well 
the upper critical field. Even for the simple model which is considered in this 
paper an extensive numerical calculation is required to obtain the details of 
the mass enhancement because it depends on the direction of the electron 
wave number k in both metals. 
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A Appendix 



A.l Connection with Green function self energy 

In the derivation of the mass enhancement in a pure bulk metal S one gen- 
erally starts from the fully occupied free electron Fermi sphere. A standard 
treatment uses the Green-function method. Here one calculates the self- 
energy of an additional electron k just above the Fermi energy (see for ex- 
ample j2], CP)- This has two contributions. Fig.4 shows the well-known 
processes involved. 

(a) The inserted electron in state k emits a phonon (q, A) and makes a vir- 
tual transition into the state k". 

(b) An electron in the state k blocks all electron-phonon processes in which 
an electron from an occupied state k' emits a phonon (q, A) and makes a 
transition into the state ko. 





Fig.4: The two contributions to the self energy and 
the mass enhancement 

The corresponding self-energies of the state k" are 



S a (k ,£) = £ k », q „ 
E 6 (ko, E) = -£ k , . 



|<7k"-ko,-q,A| 



E — Ek» — fku + iS 

I |2 
Iffko-k'-q.-M 

q,A E - £ k / + ftw + i8 
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The inverse dressed Green function has then the form 

G- 1 = E- e ko +i5-^ 

where £ = S a + £ 6 . The new quasi-particle energy is given by the pole of G, 
i.e. by solving the implicit equation 

£ k0 =£ k0 + £(k ,£ k0 ) (8) 

ii 2 ii 2 

Ifi'k"— ko,— q,A | ^ Ifl'ko— k',— q,A | 

The self-energy depends on the momentum k only through the matrix 
element, and this dependence is very weak and can be neglected. If k lies 
directly above the Fermi energy, i.e. k = k^ then we have the relation 



E k+ = E 

F 



( k ^k+) 

Next we can expand £ (k , E) in terms of E about the energy E k + 

E (ko, E) = E (k+, B k} ) + (E - B k} ) A E (k+ E k ,) 

= E (kf , B k j) - £ k + ^E (kf , £ k +) + £^E (kf , £ kf ) 
This yields for for the quasi-particle energy 

£ ko = £k0 + £ (k+ £ k +) + (£ k0 - £ k +) — Re £ (k+, £ k +) 

^o + s(4^ k ,)-£ k ^s (k+ E k .) 



-^k 



(l-^Re£(k+,£ k + 



The terms — Re £ ^kj, _E k + j +-^ k +^ Re £ ^k^, i? k + j yield essentially a con- 
stant energy shift which will be absorbed in the chemical potential. Then 
the Green function takes the form 

G(ko,E) 



E-e^- E^ Re £ (k+, £ k +) + *r ko 
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or 

G(k ,E) = 1 



z E- ^ 

where 



Z : [ l-^S!k:.l,. ] 1 > 1 



and Tk is an imaginary contribution from the self energy which we neglected 
in the discussion. 

When one performs the derivative of the self energy one realizes that its 
real part represents the realative occupation of the states k". As an example 
we take S a and obtain for — <9Re£ a (kp,E k +^j jdE. 

a I |2 I |2 

__^_v^ |gk"-k u ,-q,A| \9v."-V a -<\,\\ 

arnZ^k",q,A T? _ _ fc,J E k F - Z^k",q,A , x 2 

^ ^ £fc ^ f £ k + - £fc» - fojj 

For each k", q the right sum represents the occupation of the electron state 
k" and a phonon q due to a virtual transition from a state kp just above the 
Fermi energy. The total sum represents the reduction of the occupation of 
the state kp (before normalization). 



A. 2 Satellite states in the ground state 

At T = in the absence of the electron-phonon interaction all states within 
the Fermi sphere with k < kp (kp is the Fermi wave number) are occupied 
and all other states are empty. We denote this ground state of the electron 
system as |^ ) = II c£, |$ )- We have virtual electron-phonon emissions 

k'<k F 

from the occupied state k' into the empty states k" emitting phonons (q, A). 
Now the state k' makes a transition into one (or several) states k" and creates 
a phonon (q, A) . The resulting state can be described as 

4' = (4' + Ek'',q,A a k'',k',q,ACk'' a q,A) 

= i 1 + Ek",q,A a k",k',q,ACk" C k'a* iA ) <£, 

The amplitude of the satellites we denote as ak",k',q- This state is not nor- 
malized. 
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For the new ground state we make the product ansatz 



*o (*) 



II (l + Ek'',q,A a; k",k',q,ACk" C k'<A) l*0) 



-fa* 



(9) 



where E Q is the new ground state energy. All products over k 7 and summa- 
tions k" are restricted to k' < kp, k" > kp. The hamiltonian is 

H = Ep £ P C P C P + Ep 1 ,p a ,q'^p a -Pi,q',Ac£ a Cp 1 K',A + O* q/, A ) 

The Schroedinger equation is 

h¥ = e ¥ 

In order for ty in equ. to be an approximate eigenstate of the hami- 
tonian to first order in the electron-phonon interaction gp 2 - Pl ,q',A the states 

(°k' + Ek" q A a k",k',q,ACk»a qi A j l^o) must be (approximate) eigenstates of the 
hamiltonian, i.e. 

H ( C k + Ek»,q 1 A a; k« ) k',q ! ACk»<A) l $ o) = Vv (<& + Ek'',q,A«k'',k',q,AC k " a q ,a) l $ o) 



This yields 



Vk' (4' + Ek'',q,A«k»,k',q,AC k ,,a* )A ) |$ > 

r ^ + e k",q,A ( £k " + ^) CKk",k',q,A c k" c k'Q qi A 

= Ek' +Ek",q,A5'k"-k' -q,AC k " C k' a q,A 

. +Ek",q,Afi , k'-k",q,A c k' C k" a q,Aak",k',q,AC k //C k /a qA 

* (4' + Ek",q,A«k",k',q,AC k »a qiA ) |$ > 



This yields 



Vk' — £ k' + Ek",q,A5 , k'-k",q,A Q; k",k',q,A 
^ k 'ak",k',q,A = (^k" + faj) ak",k',q,A + <?k"-k',-q,.> 



It follows that 



«k",k',q,A 



gk"-k',-q,A 
(»& _ £ k" - hjj) 



(10) 



with the self-consistency condition 



%' - £ k' + Ek",q,A73 



<?k'-k",q,A 



(?7 k ' ~ £ k" - M 



(12) 
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There are two approximations involved in the product ansatz: (a) The 
Pauli principle excludes for c*,c*, = ^+^ tqx ^KK,^{ a ^,\) * { c k + 
the double occupancy of the state c£„c£„ for k x = k 2 . (b) The electron- 
phonon interaction can introduce a transition from the state c£ 2 into the 
satellite of c£, yielding a state c£//a* l A c£/ a* 2)A . This state is neglected. 

The total ground-state energy in this approximation is 

Eo = 

The occupation of the states c k „a* A is given by 

I |2 
, ,2 |<7k"-k',-q,A| 
|ttk»,k',q,A| - — —2 



9k"-k',-q,A 



(^ k /-£k"-^) 



1/2 



To normalize each state c£, it has to be divided by 

It should be emphasized that the energy 77^, is not the energy of a ho 
at k'. It is here only a mathematical abbreviation. 

A. 3 Ground state plus one electron 

Now we perform self-sonsistent perturbation calculation starting with the 
unperturbed ground state plus one electron k . We call this state ^o-k 
(with |k'| < k F ) 

*o;k = 4 n4' |$o) 
k' 

The resulting state in the presence of electron-phonon interaction is 

k' 

with c*, = (c* k , + Ek'', q ,A«k'',k', q ,ACk»a qiA ) where k" ^ k . 

When we derive the corresponding Schroedinger equations we obtain for 
|k'| <k F 

#k"-k' ,-q,A 
ttk",k',q,A — 7 Z r 

- ^k» - M 

?7k' = £k' + Z]k", q ,A5 , k'-k",q,Aak",k', q ,A - ^k'-ko,q,AQ!ko,k',q,A 
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The amplitude av,k',q,A has a slightly different energy denominator compared 
with the ground state since r/£, is replaced by r]^. However, the difference is 
of third order in the electron-phonon matrix element and will be neglected in 

our approximation. For the state c ko = ^c£ o + ^ k „ q x Q; k",ko,q,A c k" a q,A) one 
obtains in analogy 

<7k"-k ,-q,A 
ak",k ,q,A — 7 Z-^T 

The total energy of \Po;k becomes 

^Ojko = Ek'^k' + ?7ko 

= E + £ ko + Zlk'',q,Afi'ko-k",q,A«k",ko,q,A ~ ^k'^k'-kcq.Attkck'.q.A 

The quasi-particle energy E^ of the state c ko is then 

-Eko = ^k + Sk",q,A5 , ko-k",q,A«k",ko,q,A ~ Sk'S'k'-ko.q.A^ko.k'.q.A 

which yields 

II 2 II 2 

Tp _ i v-> |5 f k // -k ,-q,A| Iffko-kV-q^ ,-. „\ 

k0 ~ k0 + ^ k '' q < A (% - ^ - M fa* - e k0 - M 

If one compares this expression for E^ Q with the Green function expression 
in equ. (jHJ) one recognizes differences in the energy denominators. In the 
Green function expression the k -energy in the denominators is given by 
Ek . There the transition from k' to ko is replaced by a transition from ko to 
k'. The physics behind this is not obvious. As a matter of fact Schrieffer uses 
in his Superconductivity book both approaches in parallel in the discussion 
of the electron-phonon self energy. 
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